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Abstract 

A new heterodyne detection method is suggested for detecting gravitational 
waves in a Michelson Interferometer. The method is based on interference 
between phase changes which are induced by a vibrating mirror with phase 
changes which are due to the gravitational waves. The advantage of using a 
second order correlation function in the present analysis is discussed. 
In the present article a special heterodyne scheme is suggested for detecting a periodic 
gravitational wave 1-4 , which has a very low frequency (relative to optical frequencies) in a 
Michelson interferometer which has two long arms (in the order of few kilometers length). 
Since a polarized gravitational wave is expected to produce very small phase differences 
between the two arms of the interferometer , detection of such small phase differences 1-19 is 
related to optical problems. In the Michelson interferometer 20 ' 21 monochromatic light from a 
laser with a certain wavelength A is split into two orthogonal beams via a beam splitter. Each 
beam traverses a path length L and is reflected by a mirror to the beam splitter. Here the 
beams are recombined and the output is detected by a photodiode. The output observed will 
depend on the relative path difference between the two arms. In the conventional approach 
, the effective path difference is expected to be changed by the gravitational wave which will 
be observed as a shift in the interference pattern. 

In a previous study 22 it has been suggested to insert in one arm of the interferometer 
a slab of crystal (e.g. KDP) in which the index of refraction will be changed periodically 
by electrooptic effect 21 ' 23 , and the phase changes due to the modulation of the optical path 
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difference will interfere with those of the gravitational waves. This idea is however not 
feasible since the insertion of the crystal in the interferometer will introduce extra noise 
terms, and also the crystal will be destroyed by the strong laser beams which are used for 
gravitational waves detection. Since the mirrors of interferometric detectors of gravitational 
waves are often suspended by the use of magnetic forces , it is relatively easier to use 
a periodic magnetic (or electrostatic) field which would lead to periodic vibrations of one 
mirror 14 ' 19 . The modulation of the path difference between the two arms of the interferometer 
due to the vibration of the mirror can lead to interference with a periodic gravitational 
wave, which can be observed by heterodyne detection methods. The advantage of using 
a heterodyne detection method follows from the fact that it is phase dependent and can 
concentrate the measurement into a very narrow bandwidth. Since the main noise sources 
in gravitational waves detectors (e.g. shot noise, thermal noise etc.) are proportional to the 
bandwidth of measurement , the narrowing of the bandwidth by the heterodyne detection 
leads , effectively , to reduction of the noise terms. 

The interference of light from the two arms of the interferometer can be described 
schematically as 



Here A$! and A<£> 2 include all the phase terms in the first and second arms of the interfer- 
ometer , respectively. () indicates averaged values over times which are very long , relative 
to the time period of the optical wave , which has a frequency u> . G is the amplitude of the 
light output from the interferometer and Go is the the light amplitude in each of the two 
arms of the interferometer. We assume that (A$ x — A$ 2 ) includes the gravitational phase 
modulation 




(1) 
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a grav COS (u g t) 



(2) 



and phase modulation by the periodic vibration of the mirror 



A® vib = a vib cos(u vib t) 



(3) 
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By substituting Eqs. (2-3) into Eq. (1) , developing cos (A$i — A$ 2 ) with the use of Bessel 
functions of the first kind 23-24 , we find that the light intensity of the laser output includes 
the heterodyne term: 

I Het oc Ji (a grav ) sin (u g t) x J x (a vib ) sin (u vib t) (4) 



where J\ is the Bessel function of the first kind. By a straight forward calculation of Eq. 
(4) we find that the output of the interferometer will include light modulated with the 
frequencies uj 9 ± u^b- F° r Ji (i grav ) one can use the approximation: 

Here AL g is the maximal path difference between the two arms of the interferometer , 
induced by the periodic gravitational wave. For a strong periodic gravitational wave , 
an order of magnitude ^f- ~ 10~ 22 can be estimated 1-19 for a path difference L. For 6 
kilometers path difference , between the two arms of the interferometer , and A = 10~ 4 cm 
we get a grav ~ 2 x 1CT 12 which is a very small number. However J\ (a v n,) can be , relatively 
, a quite large number. Following the above analysis , we find that the intensity of the light 
output from the Michelson interferometer can have the component 

Inet oc a grav Ji (a vih ) cos [(u g + u vih ) t] (6) 

In the common method of measuring periodic gravitational waves by observing a shift in 
the interference pattern , the time of measurement is restricted to half a cycle. Since the 
flux of gravitational wave energy on earth is extremely weak , the common method has the 
disadvantage of getting extremely weak effect during this time of measurement. Another 
method in which the time of measurement is very long relative to the time period of the 
gravitational wave can be made by the measurement of the Fourier transform of the second 
order correlation function. The normalized light intensity-intensity correlation function 20 ' 25 
is given by 
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where I is the light intensity and : : denotes normal ordering of the electric field operators 
, which are in the numerator of (r). The normalized light intensity-intensity correlation 
function can be obtained by using 50-50 percent beam splitter , at the interferometer output 
, and converting the interferometer light into currents i (t) and % (t + r) by two photomulti- 
pliers. By using a spectrum analyzer and performing the Fourier transform of the current 
fluctuations one gets 20 ' 25-27 : 

oo oo 

J (<i (t) i (t + r)> - (if) e'^dr = 2eG (i) + (i) 2 J (g^ (r) - l) e'^dr (8) 

—00 —00 

Here G is the amplification factor of the photomultiplier and in deriving Eq. (8) a current 
5 (t) pulse function is assumed. 20 ' 25-27 The first term on the right side of Eq. (8) which 
is known as the shot noise 20-21 ' 25-27 can be reduced relative to the signal by using high 
intensity lasers. However , other noise terms (e.g. thermal noise, seismic etc.) will affect 
the detection of the gravitational waves. Since we are treating the detection of gravitational 
waves in the frequency range 10-10 4 Hz , there is a very strict requirement on the stability 
of the laser which should be coherent during a time r which is of order the time period of 
the gravitational wave (10 -1 — 10 -4 ) , sec. This is the main difficulty in using the present 
method but with the present advanced technology this requirement can be achieved. On 
the other hand we find that in using Eq. (8) one can average all terms over very long times 
, relative to the time period of the gravitational wave. By using Eq.(8) with the heterodyne 
term of Eq. (6 ) one may observe certain peaks in the light intensity spectrum at frequencies 
w V ib±co grav , which will indicate the observation of the periodic gravitational wave. 
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